Spatiotemporal phase synchronization in a large array of convective oscillators 
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In a quasi- ID thermal convective system consisting of a large array of nonlinearly coupled oscillators, 
clustering is the way to achieve a regime of mostly antiphase synchronized oscillators. This regime 
is characterized by a spatiotemporal doubling of traveling modes. As the dynamics is explored 
beyond a spatiotemporal chaos regime (STC) with weak coupling, new interacting modes emerge 
through a supercritical bifurcation. In this new regime, the system exhibits coherent subsystems of 
antiphase synchronized oscillators, which are stationary clusters following a spatiotemporal beating 
phenomena (ZZ regime). This regime is the result of a stronger coupling. We show from a phase 
mismatch model applied to each oscillator, that these phase coherent domains undergo a global phase 
instability meanwhile the interactions between oscillators become nonlocal. For each value of the 
control parameter we find out the time- varying topology (link matrix) from the contact interactions 
between oscillators. The new characteristic spatiotemporal scales are extracted from the antiphase 
correlations at the time intervals defined by the link matrix. The interpretation of these experimental 
results contributes to widen the understanding of other complex systems exhibiting similar phase 
chaotic dynamics in 2D and 3D. 
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1. Introduction 

Complex systems, in nature, may undergo phase synchro- 
nization transitions driven by inertia, e.g. spatially ex- 
tended multicellular systems of self-sustained oscillators 
in physics, chemistry, biology and neuroscience. The de- 
parture regime is a spatiotemporal chaos regime (STC) 
developed from a weak coupling between the units of 
the ensemble. However, under certain conditions, from 
this initial regime of STC, complex systems may display 
a simplified dynamics, which is described by fewer de- 
grees of freedom, due to the effect of synchronization sus- 
tained by nonlinear it ies. For instance, in hydrodynamics, 
for highly developed turbulent flows, the interaction be- 
tween vortices produces low-dimensional coherent struc- 
tures like the large-scale vortices observed in von Karman 
flow [de la Torre & Burguete, 2007], or for weakly devel- 
oped turbulence, far above the stationary rolls spirals 
arise in Rayleigh-Benard convection [Morris et al.^ 1993], 
During the last years, the challenging field of complex 
networks has given rise to extensive reports [Boccaletti 
et al, 2006] and renowned articles [Strogatz, 2001]. Ex- 
perimental works in the field of complex networks with 
coupled nonlinear oscillators have been performed in: 
Josephson junction arrays [Fabiny & Wiesenfeld, 1991], 
chemical oscillators in the Belousov-Zhabotinsky reaction 
[Fukuda et al.^ 2005] and in the CO oxidation on plat- 
inum crystal [Falcke & Engel, 1995], extended arrays of 
nickel electrodes in sulfuric acid [Kiss et al.^ 2002], arrays 
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of semiconductor lasers [Kozyreff et al.^ 2000]. There are 
recent theoretical models performing extended multicel- 
lular systems with ring geometry [Abrams & Strogatz, 
2006; Rajesh & Sinha, 2008]. 

Depending on the range of interactions between the 
oscillating units in a network, synchronization phenom- 
ena can be classified into: global synchronization (all-to- 
all), local synchronization (between nearest neighbors) 
and nonlocal synchronization (farther than the nearest 
neighbors). There are recent research works interested 
in time-dependent phenomena driven by nonlocal cou- 
pled oscillators [Tanaka & Kuramoto, 2003; Shima & Ku- 
ramoto, 2004]. A priori, under critical controlled condi- 
tions interactions between convective diffusively coupled 
cells are expected to be local. 

The dynamics of complex system consisting of an en- 
semble of coupled oscillators, far enough from the first 
instability, can be described from their phases (j)i{t). This 
phase dynamics can shed light into the understanding of 
critical phenomena when a subsystem of oscillators bi- 
furcates towards a new oscillatory state represented by 
an average phase ^{t). This synchronized ensemble of 
oscillators defines a cluster. 

First attempts on phase synchronization or self- 
entrainment theory were originally introduced by Win- 
free for biological communities [Winfree, 1967], and af- 
terwards this theory was developed by Kuramoto [Ku- 
ramoto, 2003] in the phase reduction theory of weakly 
coupled limit cycles (KM) for an ensemble of N identical 
oscillators with global coupling. To characterize a phase 
synchronization transition, we look for an order parame- 
ter (or parameters) that quantifies the way in which the 
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system evolves towards the new state. From the aver- 
age synchronized phase of the system <l> (cluster phase 
synchronization), KM [Kuramoto, 2003] defines the or- 
der parameter re**^ for the ensemble of N oscillators with 
random phases (j)j : 

re'^ = iVe'^^' (1) 

3 

We report experimental results concerning a quasi- ID 
convective system whose behavior can be understood as 
the resulting dynamics of A/" 80 convective cells {os- 
cillating thermal plumes). These oscillators define a ge- 
ometric network as they are placed on a spatially ex- 
tended array. Coexisting with a stationary pattern (ST) 
of wavelength kg , which imposes a discrete space symme- 
try, there are irregular and stationary clusters of time- 
dependent patterns where oscillators are preferably in 
antiphase (counterphase) with their nearest neighbors, 
the kind of limit-cycle oscillators. Similar patterns of 
clusters have been observed in the Belousov-Zhabotinsky 
reaction-diffusion system with global feedback [Vanag 
et al, 2000]. 

The results reported here contribute to study time- 
varying topologies as we increase the driving force. We 
extract the topology from the dynamics, and we charac- 
terize the synchronization clustering process in the spa- 
tiotemporal beating regime (ZZ) from the antiphase ma- 
trix. Regarding topology, we find parallel contributions 
to the field of complex networks [Timme, 2006; Bialon- 
ski & Lehnertz, 2006]. Other efforts towards the under- 
standing of similar systems are developed in terms of dy- 
namical weights from adaptive processes [Zhou & Kurths, 
2006] , taking into account a varying coupling strength. 

The aim of this paper is to understand the natural phe- 
nomenon of interacting oscillators in a ID convective sys- 
tem. Inside the regime of spatiotemporal beats, the cou- 
pling between oscillators becomes nonlocal as the driving 
force is increased. This nonlocal coupling is responsible 
for a global transition deduced from the correlation times. 
For each control parameter value a time-dependent con- 
nectivity matrix is built. An order parameter is obtained 
from the cross-correlations between the critical phases 
{(j)^) of N oscillators, and it is shown to be of the order 
of the spreading size of the new synchronized phase. 

2. The Convective Problem 

Our convective layer of fiuid is situated in a rectangu- 
lar vessel {Lx x Ly, Lx ^ Ly) opened to the atmosphere. 
The fiuid used is transparent to visible light (silicone oil 
with a viscosity of 5 cSt). Quasi- ID-dynamics is caused 
by heating the fiuid layer along a line which is placed 
underneath and at the center. This line is provided by 
placing a heater rail below the mirror that sustains the 
fiuid layer. The temperature at the lateral walls in the 
larger direction is controlled at 20.0 dz 0.1 °C, the same 
one as the room temperature. The array of convective 
cells is observed using the shadowgraphy technique which 



reveals the fiuid dynamics inside the bulk of the layer on 
a screen. The control parameter is the vertical temper- 
ature difference, AT, between the constant temperature 
along the heating line and the room temperature. Ex- 
tended details on the experimental setup are found in a 
previous work [Miranda & Burguete, 2008]. We have ex- 
plored the range AT = [0 — 30] for depths of the fiuid 
layer of d = [2.5 — 5] mm. The corresponding stability 
diagram {d vs AT) has been already reported from the 
analysis of the spatiotemporal signals S{x^ t) [Miranda 
& Burguete, 2008; Miranda & Burguete, 2009]. At each 
asymptotic state, these signals are obtained by recording 
the shadowgraphy image along a line x = [0, Lx] next to 
the center. 

In the shadowgraphy image [Fig. [ija)], an oscillator 
will be associated to each one of these convective cells 
(or hotspots) along the heating line. The acquisition line 
to record the spatiotemporal signal S{x^t) is represented 
in Fig. [TJa) by a blue line. This line is placed parallel 
to the aligned uprising hotspots in the center. The cor- 
respondent brightness profile [Fig. mb)] allows to track 
the oscillators position in time. This is given by selecting 
the implied relative extrema from Fig. [Tfb). We choose 
a characteristic spatial scale defined by the wavelength 
2A5 of the oscillating thermal cells in these instabilities. 

In Fig. [2ja) we show a typical spatiotemporal diagram 
S{x^t) above the threshold of the secondary bifurcation 
to a zig-zag regime, the ZZ pattern. This spatiotemporal 
diagram corresponds to the acquisition of global informa- 
tion about the interaction of N elements upon time. This 
global information, as it is explained later on in Sec. 3, 
can also be divided into N-oscillating units [see the corre- 
sponding tracking of N = 5 oscillators in Fig.[2l^b)]. From 
each S{x^t) diagram, bidimensional Fourier fast trans- 
form captures the essence of the global oscillations into 
a set of regional peaks. In the Fourier space, each peak 
corresponds to a mode with a well-defined "spatiotem- 
poral frequency band" . From the data analysis, for each 
fundamental mode we get M(A:,cj), where the wavenum- 
ber k and the frequency 00 are the maxima of each peak 
[Miranda & Burguete, 2008; Miranda & Burguete, 2009]. 

For a constant depth d = 3 mm, as AT is increased 
beyond the ST regime, waiting for the asymptotic state 
to settle, the system goes through a "mixed pattern" in 
a STC regime. This STC regime is characterized by the 
presence of several irregular clusters. Oscillators belong- 
ing to this clusters perform an alternating pattern (ALT) 
[with wavenumber ks/2 = ks/2 and with frequency uOait] 
which coexists with the ST pattern in a mixed ST/ ALT 
regime. From the coalescence of these various spreading 
phase subsystems the dynamics undergoes a supercritical 
bifurcation towards a new state of spatiotemporal beats 
(ZZ). This ZZ pattern is composed of two phase sub- 
systems with high spatiotemporal coherence. It is char- 
acterized by a zig-zag geometry due to the splitting of 
the temporal frequency into (oOait^^zz) with cOzz ^ ^ait- 
Also, there is a spatial frequency splitting into (A:s/2, kzz) 
with ks/2 ^ kzz, which is responsible for the existence of 
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two high coherent patterns. This ZZ regime cahs for a 
phase description obtained from the signal of each os- 
cillator. At the STC regime, weak coupling can sustain 
irregular clusters meanwhile, in the ZZ regime, a stronger 
coupling locks the phase into stationary clusters. From 
the diagrams we have shown in a previous work 

[Miranda & Burguete, 2009] the evolution of the funda- 
mental modes in the mixed ST/ ALT and the ZZ regimes, 
and also that the oscillating frequencies increase linearly 
with the control parameter. For thicker layers the oscil- 
lators are able to synchronize moving towards one privi- 
leged direction, a travelling wave regime (TW) [Miranda 
& Burguete, 2008]. 

Phase mismatches between 50 oscillators, calculated 
using a spatiotemporal cross-correlation signal process- 
ing, provide the framework to understand the spatiotem- 
poral coherence of the beating regime from a nonlocal 
coupling. New dynamical conditions like stronger non- 
linear coupling set a longer interaction range, further 
than first neighbors (three interacting oscillators). An 
increased number of interacting oscillators is the result 
of strongly interacting convective cells, we might think of 
a new structured flow that encloses several long-ranged 
interacting cells. 

On the other hand, connection topologies are not fixed 
because the topology of the network is defined by time- 
varying links between the oscillators. The loci of the links 
are determined from the collapsing interactions between 
the nearest neighboring oscillators (link matrix). In these 
contact interactions two convective cells are colliding. 

Also, we measure the synchronization range (number 
of synchronized oscillators) given by the correlated syn- 
chronized phases These phases are obtained demod- 
ulating the raw phase of each oscillator nearby the 
selected synchronization frequency {coi = dt(l)i). These 
raw phases are functions of the several unstable or 
critical temporal modes. 

Hence, in our ID-array, a "critical phase" (^^) emerges 
at the threshold of the synchronization transition for each 
oscillator. This critical phase (j)^ slightly shifts from the 
cluster phase ^. Under these conditions, the developed 
analysis techniques reported here are somehow inspired 
by the Eq.( [T]) (KM). The key is to build a minimum 
cross-correlation matrix (antiphase matrix) and the cor- 
responding matrix of time lags at minimum correlation 
(= maximum antiphase correlations) between each pair 
of oscillators. 

3. Spatiotemporal Phase Synchronization Tran- 
sition Scenario 

In our ensemble of oscillators, the phase chaotic regime 
extends over the STC and ZZ regimes. Firstly, short- 
range interactions (due to diffusive coupling) are able to 
develop a regime of STC. In this regime, spatiotemporal 
correlations in the stationary pattern decay due to the 
increasing presence of irregular clusters. This chaotic dy- 
namics corresponds to an inhomogeneous growth of the 
ALT pattern until a critical value [Miranda & Burguete, 



2009]. Secondly, a phase synchronization regime is char- 
acterized by the presence of beats, ZZ (see the sketch for 
the temporal beats in Fig. EJc)). These beats are due 
to a small shift in the spatial and temporal frequencies. 
From now on, we will focus on the latter synchronization 
transition. Experimental results follow this ZZ regime for 
an ascending sequence of 5 K by steps of 0.5 K. Spatial 
desynchronization appears for higher control parameter 
values with the presence of phase singularities. Thus, a 
temporal synchronization only remains. 

When the original signal along the entire ensemble 
S{x^t) is decoded for N = 50 oscillators [e.g. for a 
few oscillators from the N-ensemble see Fig. Efb)] we ob- 
tain an individual signal for each i-oscillator: A'i(t) = 
Xi + Ai cos[^ i{(j)f^ ^ (j)^^ ^ (j)^)]^ where is the raw phase 
containing all the involved critical phases (j)^ (for the ALT 
and ZZ characteristic phases), and 0° is a constant phase 
shift. In this way, we are able to understand how do the 
oscillators interact in space and time. 

We demodulate around the emerging unstable mode 
cof^ = dtcj)^^^ the new phase associated to each oscillator 
is simply given by 4'i{= (t^l^)- The cross-correlation phase 
vector is given by: 

C.,(c^..,rO = (cos((/)f )cos((/)f )) (2) 

where ti are the time lags, e.g. if we take a couple of os- 
cillators (i, j) in antiphase, this behavior belongs to the 
characteristic alternating pattern ALT (two counteroscil- 
lating modes with icjazt), then ti is given by the half of a 
period T (T/2 is the time that takes the z-oscillator Xi(t) 
to catch up the j-oscillator Xj{t + T/2)). Because the 
great majority of oscillators are in antiphase, the cross- 
correlation vector Cij happens to show the strongest cor- 
relation between neighboring oscillators in the minima of 
the cross-correlation vector (antiphase conditions). This 
fact naturally comes from the type of attractive-repulsive 
interaction between oscillating units. In order to extract 
information about the phase correlations along the array 
we study this antiphase matrix [see Fig. [3l^a,b)] defined 
as: 

Aij{uJzz) = min Cij{uJzz,ri) (3) 
n 

This is a TV X matrix whose elements are given as 
the minimum of Cij (maximum antiphase correlations) 
at their corresponding time shifts r^ Conclusively, we 
may consider that this matrix is actually supplying the 
"degree of phase shifts" between each oscillator regard- 
ing the cluster synchronization phase: — <l>|. In this 
way, the cluster phase ^ fulfills the following definition 
dt^{t) = cOzz = XliLi^f^/^' where uOzz is the average 
critical mode of n < synchronized oscillators. As 
we show in Fig. [3Ka,b) the surface representation of the 
matrix Aij{ujzz) contains the coupling nature between 
all-to-all oscillators of a ID array. In these figures the 
depressed regions (the deeper ones) correspond to the 
higher synchronized oscillators defining the boundaries 
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that surround clusters. Hence as we increase the control 
parameter a new phase synchronization domain emerges 
which correspond to the spatial beating phenomena: two 
stationary clusters [Fig. [3jb)]. 

Along the ascending sequence in AT, from the analysis 
of Ti applied to an oscillator belonging to the cluster (n = 
31), it has not been possible to establish any dependence 
on AT. From this point on, our research follows the 
information given the antiphase matrix, Aij((jJzz) whose 
values inside a cluster keep sufficiently small, below a 
fixed correlation value Aij{ijOzz) = e (for this sequence 
e = —200 [Fig. [3l(a,b)]). Therefore, e is kept constant 
in order to guarantee the phase locking between neigh- 
boring oscillators for the whole oscillators array and for 
a whole sequence of measurements. When the first sta- 
tionary cluster emerges from AijiuOzz) the critical control 
parameter is AT^ = 31 K [Fig.IS^a)]. 

We define the order parameter Z^uOzz) to quantify the 
dynamics of the complex network as we increase the con- 
trol parameter. Firstly we define for each z-oscillator 
^^ii'^zz) (as we move from the left to the right): 



The result given by Ziiujzz^ = ^ means that for a given 
i-oscillator, the following number rj < N of oscillators 
are highly correlated to it. From this analysis we are 
able to determine the maximum number of synchronized 
oscillators in a cluster at the asymptotic state determined 
by AT. Thus, the order parameter is obtained from: 



Z{uJzz) = max {Zi {ujzz)) (5) 

for an z-oscillator belonging to the cluster. In Fig. [H 
we find out, far from the threshold (AT > 33 K) and for 
the spatial coherent subsystem (stationary cluster in ZZ), 
that 15-17 oscillators are coupled, but for the same pat- 
tern and lower control parameter values approximately 5 
oscillators are coupled. 

The outstanding presence of the spatial coherent do- 
mains (spatial beats) as the system moves forward a 
global synchronization instability is shown from the av- 
erage correlation time (r) [Fig.[5](a)]. The assigned aver- 
age value (r) for each control parameter value is obtained 
from the correlation times of the whole array r^. These 
are measured from the maximum time interval at which 
\Cij{uJzz,ri)\ decays to \Cij{uJzz,ri)\e-^ for ri = 0. In 
Fig.[5fb) we show the distribution of for three represen- 
tative control parameter values. Far from the threshold 
(at AT = 34.5 K), reaches the top time recording value 
for almost all the oscillators, in contrast to the values ob- 
tained below the threshold (at AT = 28 K). Clearly, at 
the threshold (ATc) the oscillators belonging to the clus- 
ter show the longest correlation times [Fig. [5l^b)]. 



But there is still one missing piece in this network puz- 
zle, when we ask about what is going on with the inter- 
action range in the overall ascending sequence towards 
this global synchronization instability. Regarding a non- 
linear hydrodynamics point of view, we might think of a 
new dynamics of the convective cells in the phase chaotic 
regime, a kind of package An consisting of more than 
three interacting oscillators. This wider "enveloping con- 
vective cell " could not be the effect of a merely diffusive 
local coupling as for the ST regime. The ZZ regime might 
be sustained by a stronger nonlinear coupling between os- 
cillators. If we choose three oscillators from the array [see 
Fig. EJa-c)]: n = 25 (this oscillator belongs to the front 
of the right hand side cluster only far from the thresh- 
old [Fig. [6jc)]) and n = 35,37 (these ones belong to the 
cluster on the right hand side [Fig. [6l^b,c)]), then repre- 
senting \Anj{uJzz)\ for a fixed i = n and j = 1,...,A^, 
from Fig. El^a-c), we find out higher phase correlation 
values, and therefore stronger coupling, the further the 
control parameter is from the threshold. Besides, the 
decay tendency of the antiphase cross-correlations be- 
tween the selected oscillators and the nearest neighbors is 
smoother far from the threshold [Fig. [6l^c)], considering 
the abrupt decay of \Anj{uJzz)\ between nearest neigh- 
bors in Fig.ini^a). We observe how this strong coupling is 
already outstanding at the threshold [Fig. [6^b)]. These 
last results might allow us to determine an effective non- 
local interaction range An from the cross-correlations. 
Nevertheless, the topology (link matrix) will be involved 
in these results, as it will be shown in the next section, 
from the colliding positions between oscillators. 

4. Time- Varying Topology 

Our objective is to decode some topological aspects such 
as geometric loci of the oscillators and the degree of con- 
nectivity from the dynamics of our ID coupled network. 
It is interesting to consider the degree of attractive inter- 
action between oscillators, regarding this fact we could 
define two kinds of connectivities [sketch in Fig. [2jc)]: 
type I for an overlapped pair of oscillators, a collision 
between two adjacent hotspots has been produced; and 
type II when the distance between two neighboring oscil- 
lators reaches a minimum but without collision. In the 
following, connectivity always refers to type I. 

In the synchronized regime ZZ, we build the link ma- 
trix Lij by locating in the tracking image [e.g. Fig. EJ^b)] 
the contact loci between two adjacent antiphase oscilla- 
tors that satisfy: A'^_i(t) = A'i(t) [see sketch in Fig.[2fc)]. 
A link point is a collision position, this means that a cou- 
ple of oscillators overlap at the same time that the linking 
is affecting at least 4 oscillators. After filtering contact 
links provided only by oscillators with frequency cOzz^ Lij 
is nonzero only at these contact loci [e.g. Fig. [3jc,d)]. 
The analysis of L^j shows that as the system is getting far 
from the threshold the number of links increases. More- 
over, it allows to study how the dynamics of the rewiring 
connections, for different control parameter values AT, 
plays a role in the coupling interaction between oscilla- 
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tors. 

The link matrix Lij provides quantitative information 
about the evolution of the number of links in time. This 
information is simplified by taking running averages over 
a time interval [Fig. [7]. We choose 62 s as the time to 
average, which approximately corresponds to half of the 
periodicity of the temporal beating. From Lij along the 
ascending sequence, we recover the characteristic time 
periods from these graphics [Fig. [7]. The envelope of the 
temporal beat and the internal wave (carrier) periods are 
approximately 120 s and 15 s respectively [Fig. [71(a,b)]. 
At this point, we look for an homologous cross-correlation 
phase behavior at times belonging to the minima ( tm 
with a lower degree of links), and maxima {Im with 
a higher degree of links) of the running averages over 
Lij. In order to study this effect we obtain the cross- 
correlation phase vector for a fixed time: 

Cij{uJzz,Xl) = {cOs{(j)i{tm,M)),COs{(l)j{tm,M))) (6) 

where xi is the oscillators lag xi = 1, • • • , A^. When the 
number of links is zero the "spatial" cross-correlation 
phase vector is chaotic, meanwhile when the number 
of links is nonzero the typical cross-correlation vec- 
tor is similar to Fig. [8](a). From the absolute value 
\Cij{u;zz^ Xi)\tm we identify a correlation length of 4-5 os- 
cillators [Fig.[8](b)] with low peaks at the cluster size (13- 
17 oscillators belonging to the clusters). Otherwise, from 
the absolute value \Cij{uJzz^Xi)\tM we identify the maxi- 
mum correlation "length" at the cluster size of approx- 
imately 15-18 oscillators [Fig. [81(c)]. But certain cross- 
correlation behavior, like the one showed in Fig. [Hd), 
still shows the maximum coupling strength around 4 os- 
cillators. Thus, far from the STC regime, in the ZZ phase 
synchronization transition, we are allowed to consider 
that the dynamics affects the underlying topology of the 
network. 

From these results, it follows that nonlocal coupling 
existing along the ZZ regime could be possible. This cou- 
pling might go beyond first neighbors (An > 3), differing 
from the kind of coupling that it was already shown below 
the threshold in Fig. [BJa) for An = 3. In the ZZ regime, 
because the interaction range is wider than a purely dif- 
fusive one we might certainly expect roughly An = 4 — 5 
oscillators. 

5. Discussion and Conclusions 

We have reported the study of a synchronization tran- 
sition beyond spatiotemporal chaos in a hydrodynamic 
system. This consists of an ensemble of N ^ 50 identical 
oscillators (the size of each oscillator is approximately of 
6 mm in a real array 80 oscillators), developed in a fluid 
layer under ID thermal convection [Miranda & Burguete, 
2009]. These coupled z-oscillators {i = 1,...,7V) whose 
signals are given by: Xi{t) = Xi+Ai cos[^i((/)f ^ (/)f^, (/)°)], 
are placed in the array at their respective "zero positions" 
(xi), which settle two very close wavenumbers: ks/2 and 
kzz] and have two close critical frequencies: cOait and 



^zz = dt^^ respectively. We have described the syn- 
chronization process in terms of the individual phases 0^ 
regarding the cluster phase that is, the analysis of the 
antiphase matrix Aij{u;zz)- The surface of AijiyUOzz) de- 
fines the stationary cluster domains for the critical mode 
ujzz-, in addition it shows no significant correlation for 
uOait' The number of synchronized oscillators is defined 
by ZiiuOzz) = T] < N for i = 1, . . . , N and agrees with the 
assigned high time correlation values r^. 

We might notice that from the time lag analysis be- 
tween each pair of oscillators (i,j), defined as ri^i = 
Tij {mod Tzz) (module Tzz = 27r /ujzz)^ is more difficult to 
recover the number of synchronized oscillators {Zi{u;zz)) 
than from the analysis of the antiphase matrix. Further- 
more, the averaged value of time lags (in the ascending 
sequence) between an oscillator belonging to the cluster 
regarding the remaining oscillators, show no connection 
with the increasing time correlations. 

The values of Zi{LJzz) for an i-oscillator belonging to 
the cluster on the left are less well-defined than the ones 
from the cluster on the right, and the corresponding am- 
plitude inside the cluster happens to be "turbulent" for 
the critical mode cOzz- It should be taken into account 
that the order parameter Z{uJzz) depends on the phase 
distribution along the ID-array. From Fig. [4] we show 
that the number of coupled oscillators belonging to the 
cluster increases from 5 to 17 for higher values of the 
control parameter. Although the frequency is a mono- 
tonically increasing function with temperature uOzziT) 
[Miranda & Burguete, 2009], Z{uJzz) is not a continuous 
function of AT . At AT = 33 K, the order parameter and 
average time correlation correspond to a desynchronized 
pattern. 

Results on Fig. EJa-c), obtained from the antiphase 
matrix, shows the type of the coupling interaction be- 
tween the oscillators at the boundary and the oscillators 
belonging to the synchronized domain. In Fig.[6l^b), it is 
shown that for the most distant oscillator (n = 25) from 
the boundary of the synchronized ZZ pattern, the phase 
correlation decays "more quickly" , as it is expected from 
a diffusive coupling of the type Anj{iJzz) ^ e~^^ (short 
interaction range). Meanwhile for the remaining oscilla- 
tors belonging to the cluster, the antiphase correlation 
has a smoother decay. This fact should be interpreted 
as the consequence of a longer interaction range in an 
"envelope roll" of An = 4 — 5 interacting oscillators. 

In our spatially extended array, we have shown the 
dependence of the network topology on the dynamics 
[Fig. [7]. The increasing number of contact interactions 
provide a sensitive dynamics which is revealed by the 
peaks in the spatial correlations between oscillators at 
\C^j{^zz,Xl)\trr^,M [^Ig- [D^a-d)] , also we get the size of the 
clusters. 

We have shown that the wiring of connections from 
the link matrix T^j, is not fixed, neither in space nor in 
time. From the phase description theory (KM), a suitable 
model for an ensemble of identical z-oscillators, in the 
framework of nonlocal two-way coupling for attractive- 
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repulsive interactions, could be (i = 1, . . . , A^): 

N 

ci>, = F,((/>f , . . . , (/>°) + ^ H,,{t) • r,,((/>, - (7) 

Each z-oscillator has its own dynamics represented 
by Fi as a function of the critical phases. The time- 
dependent connectivity function is the adjacency matrix 
^iji'^) — ^ij ' ^iji^)' From An (given by the analy- 
sis of the running averages of the link matrix Lij) we 
can build an adjacency matrix at zero time M^j(O) such 
that Mij = 1 if |i - j| < An, and Mij = other- 
wise. The time-dependent connectivity, extracted from 
Lij, might be mapped into Gij{t). The coupling func- 
tion between oscillators is Tij{(j)i — for this type of 
model equation Eq.([7j) comes from a slow perturbation 
(e.g. with Tij{(l)i — (pj) — sm{(j)i — (/>j) [Kuramoto, 1975]), 
the real coupling function in the ZZ regime might have 
a more complex behavior according to a varying cou- 
pling strength that depends on the control parameter 

r,, ((/),-(/),, AT). 

For a supercritical bifurcation, a phase description 
might not necessarily hold for any regime far from the 
first instability, because the possible coupling between 
the amplitude and the phase should play a role. Regard- 
ing this fact, synchronous behavior in the array is re- 
lated to the existence of an unstable traveling mode ujzz 
with enough amplitude to start synchronization (keeping 
strong values in \Aij{u;zz)\) sustained by nonlinearities. 
When this critical mode cOzz is missing the rewiring of 
connections is only due to the mode uJait- This mode rep- 
resents weak coupling for counteroscillating modes, and 
moreover the corresponding interaction range still rep- 
resents the kind of diffusive transport (of heat and vis- 
cous momentum) that one may expect from a local cou- 
pling [Fig. [61(a)]. In conclusion, the reported results for 
an spatially extended system show how the amplitude of 
the unstable mode plays a role in the phase synchroniza- 
tion process (e.g. according to numerical works [Vadi- 
vasova et a/., 2001]). Besides, nonlocal coupling could 
explain similar transitions in nature and experiments, as 
it has been reported in recent theoretical works [Abrams 
& Strogatz, 2006]. This synchronization phenomena is a 
new challenge that might need to be considered in the 
theory of cluster formation. 
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FIG. 1: (a) Photograph of the shadowgraphy image, dark 
spots at the center hne correspond to the ascending hot 
plumes (hotspots). The continuous blue line represents 
the acquisition line recording the spatiotemporal diagrams 
S{x,t). (b) Instantaneous brightness profile along the ID- 
array (over the acquisition line). As is the wavelength of the 
stationary ST pattern. The bright spots (resp. dark spots) 
in the top (resp. bottom) figure correspond to the position of 
two adjacent oscillators. 




FIG. 2: (a) Shadowgraphy image at AT = 32 K; (b) tra- 
jectories of the oscillators tracked from the frame in (a); (c) 
sketch of the trajectories for a few oscillators separated at 
the characteristic spatial scale 2As in the ZZ pattern. Black 
filled circles indicate the contact interaction loci (connectiv- 
ity type I) , white filled circles indicate the minimum distance 
between oscillators without collision (connectivity type II). 
Discontinuous lines are a guide to the eye showing the ob- 
served characteristic zig-zag (ZZ) pattern (temporal beats). 
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FIG. 3: Surface views of the antiphase correlation matrices 
Aij((jJzz) at the critical frequency mode Uzz for two phase 
synchronization states at (a) AT = ATc = 31 K and (b) 
AT = 34 K. The corresponding link matrices Lij generated by 
contact interactions between oscillators are shown in pictures 
(c) for AT = ATc = 31 K and (d) AT = 34. 
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FIG. 4: Order parameter values Z{uJzz) for the ascending 
sequence. The vertical dashed line split two different regimes. 
Horizontal solid lines show averaged values of the number of 
oscillators involved in clusters. 
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FIG. 5: (a) Average correlation time (r) vs AT, discontinuous 
line is a guide to the eye; (b) correlation time (n) for each 
oscillator at AT =28, 31, 34.5 K. Top and bottom continuous 
lines represent the maximum and minimum time correlations, 
dotted lines for AT =31 K are a guide to the eye. 
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FIG. 6: Antiphase cross-correlations Anj{oOzz) (j = 1, . . . , A^) 
for oscillators n = 25, 35, 37 at: (a) AT = 28 K; (b) AT = 
31 K and (c) AT = 34.5 K . Continuous and dashed lines are 
a guide to the eyes. Vertical dashed lines define the cluster 
domain. 
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FIG. 7: (a) Running averages over 62 seconds of the link 
matrix Lij at AT =31,32 K; (b) running averages over 62 
seconds of the hnk matrix Lij at AT =34, 34.5 K. Red and 
blue doubled arrows show the temporal beat period (envelope) 
and the internal wave (carrier) period respectively. 





FIG. 8: (a) Cross-correlation vector at time of maximum 
linking, Cij{oOzz,Xi)\t=404s (t — 404 s time of maximum 
linking determined by the maxima of link matrix Lij) for 
AT = 34.5 K . Module of the spatial correlations: (b) at 
minima \Cij{uJzz,Xi)\t^iQQs for AT = 33.5 K; (c) at max- 
ima \Cij{uJzz^Xi)\t^797s for AT = 34.5 K; (d) at maxima 
\Cij{uJzz,Xi)\t^Qis for AT = 33.5 K . 



